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Table 1 


6'°C and 6'N values for vent and non-vent organisms and 
their potential food sources. Sample size (n) for animals 
= the number of individuals analyzed. Three samples of 
vent bacteria and one algal sample were analyzed. Values 
are given as the mean + SE if n is greater than 1. 


Organism n 6°C ôN 

Vent 

Bacteria 3 =i (0) se 0.3} —4.1 + 1.3 

Lottia limatula 27 — 30.4 —1.3 

Mytilus edulis 3 —16.6 + 0.7 +6.4 + 1.0 
Non-vent 

Benthic algae 1 S51 +6.8 

Lottia limatula 3* = 155 +6.3 

Mytilus edulis 1 z153 +8.6 


* These limpets were analyzed as one sample. 


the nearest detectable vent opening) sites at White Point 
in southern California. Vent openings were easily detected 
because the rocky substrate surrounding an opening is 
devoid of algal growth but covered with dense, off-white 
mats of sulfur bacteria. 

Vent bacterial and non-vent benthic algal samples were 
scraped from rocks. Limpets were collected and kept live 
in 25-um filtered seawater to allow for gut evacuation, so 
the isotopic composition of the whole body could be de- 
termined without contamination from ingested but unas- 
similated material. Mussels were collected and their foot 
and mantle tissues were dissected out for analysis, to avoid 
contamination by gut contents. All samples were lyophi- 
lized, ground to powder, treated with 1 N HCI to remove 
carbonate, and lyophilized again before combustion 
(NoRTHFELT eż al., 1981). The resulting CO, and N, were 
separated cryogenically and their volumes determined 
manometrically prior to mass spectrometric analysis. The 
isotopic compositions are expressed in the usual 6 notation 
where 6X = [{R(sample)/R (standard)} — 1] x 1000%o 
and X = C or 5N, R = ®C/"C or 5N/*N, and the 
standard = PDB belemnite or atmospheric N,, respec- 
tively, for carbon and nitrogen (Rau, 1985). 


RESULTS anD DISCUSSION 


The carbon isotope ratios reported in Table 1 suggest that 
the diet of intertidal vent limpets consisted predominantly 
of sulfide-oxidizing bacteria. The 6°C values of vent lim- 
pets match those of vent bacteria. The 6°C values of non- 
vent limpets, in contrast, were decidedly less negative than 
those of the vent limpets, and are similar to those of their 
likely food source, photosynthetic algae encrusting the rocks 
on which they live. 

Carbon isotope ratios of mussels living near and away 
from the intertidal vents (Table 1) suggest that the vent 
mussels are using only small amounts of vent bacteria as 


a food source. The 61°C values for vent mussels are only 
about 1%o more negative than those for non-vent mussels. 
This 1%o difference contrasts with the 16%o difference in 
the 6'°C values between the vent bacteria and the non-vent 
algae, which may be taken as the difference between the 
chemosynthetically and photosynthetically derived food 
sources in the area. If the photosynthetically derived food 
used by the mussels had 6'°C values similar to the benthic 
algae we analyzed (Table 1), the mussel 6'°C values suggest 
that more than 95% of the vent mussel food was not derived 
from vent bacteria. Water column concentrations of sus- 
pended bacteria detached from vent bacterial mats were 
not measured in this study, but are presumably high near 
vents because many suspended filaments of mat bacteria 
are visible to the naked eye over vent openings (especially 
when wave action detaches parts of the attached mat). 
Presumably, a suspension feeder, like Mytilus edulis, that 
is situated directly over a vent must entrain some bacteria 
in its feeding current. Most of these bacteria may subse- 
quently be sorted out and rejected in favor of plankton in 
non-vent water that mixes with vent effluent. 

Nitrogen isotope measurements may also provide in- 
formation on trophic relationships. For deep-sea vent and 
other marine organisms, the ôN values increase approx- 
imately 2-3%o with each change in trophic level (RAU, 
1985). Our measurements for intertidal vent limpets are 
in good agreement with this trend. The 6'5N value for 
limpets is 2.8%o higher than that for vent bacteria (Table 
1), strengthening the conclusion that the bulk of vent limpet 
diet consists of vent bacteria. The unusually low 6°N value 
for vent bacteria may be an indication of bacterial N, 
fixation as a vent nitrogen source, rather than a biologically 
cycled nitrogen source (e.g., nitrate), which would have a 
more positive 64°N value (RAU, 1985). On the other hand, 
it is also possible that vent water contains an !°N-depleted 
organic source of nitrogen. For non-vent control limpets, 
the ôN values are close to that of their algal food source 
but do not show the slight rise predicted with a change in 
trophic level. Nevertheless, the values are consistent with 
a non-vent photosynthetically generated food source (RAU, 
1985). 

The ôN values for vent and non-vent mussels are con- 
sistent with our conclusion that vent mussels use vent bac- 
teria as only a minor food source. 

According to STEIN (1984), mollusks are the only ma- 
croinvertebrates occurring commonly at subtidal White 
Point hydrothermal vents. At the littoral zone vents ex- 
amined in our study, mollusks are the dominant macroin- 
vertebrates, but two arthropods, the striped shore crab 
Pachygrapsus crassipes Randall, 1839 (Figure 1d) and the 
blue-clawed hermit crab Pagurus samuelis (Stimpson, 1857), 
are also common. It is therefore possible that these and 
other shallow-water marine macroinvertebrates, in addi- 
tion to file limpets and black abalone, use vent bacteria as 
an important food source. Also, the assemblage of tiny 
organisms associated with mats of coastal vent bacteria 
(e.g., copepods, flatworms, polychaetes, flagellates), and 
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the bacteria themselves, have not yet been fully charac- 
terized. Because of easy access, there is great potential to 
investigate further the nature of trophic relationships with- 
in these shallow-water hydrothermal vent communities. 

Our study is a first account of both carbon and nitrogen 
stable isotope measurements indicating that geothermally 
driven production of biomass is of major nutritional im- 
portance for certain consumers, even in the sunlit zone 
where photosynthesis dominates. 
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Abstract. ‘The West Indian topshell, Cittartum pica, is the most intensively fished intertidal marine 
snail of the West Indies. However, its population dynamics are poorly known. Different habitat types, 
varying in degree of wave exposure, were chosen to study survival, growth, and fecundity of unexploited 


Bahamian Cittarum populations. 


Several aspects of Cittarium life history differed between habitats. Snails inhabiting wave-exposed 
habitat had lower rates of survival and growth and a smaller average shell size than snails of sheltered 
habitat. Wave-exposed habitat had higher Cittarium densities and higher predatory gastropod densities 
than more sheltered habitats. Dead shells recovered from wave-exposed habitat showed a higher pro- 
portion of lethal damage due to shell-drilling and crushing. Size of maturation was smaller and fecundity 
was lower for snails residing in wave-exposed habitat. 

Differences in population structure between habitats were attributed to differences in survival and 
growth and not to temporal differences in recruitment. Low Cittarium abundance in sheltered habitat, 
despite high survival rate, suggested a lower level of recruitment to such habitat. 


INTRODUCTION 


The West Indian topshell, Cittarium pica (Linnaeus, 1758) 
(=Livona pica) has been harvested since pre-Columbian 
times (GOULD, 1971) and is the largest shallow-water 
trochacean gastropod of the western Atlantic, often reach- 
ing a shell diameter of 10 cm. This herbivorous snail is 
sedentary and is found in a variety of rocky intertidal 
habitats throughout the West Indies. The snail has rarely 
been found in southern Florida (ABBOTT, 1976) and suf- 
fered recent extinction in Bermuda (CLENCH & ABBOTT 
1943). 

Information on the biology of Cittarium is limited. The 
snail’s anatomy (CLENCH & ABBOTT, 1943; GRAHAM, 
1965), shell ultrastructure (WISE & Hay, 1968a, b), and 
flight response (HOFFMAN & WELDON, 1977) have been 
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described. Only one study has dealt with Cittar1um pop- 
ulation dynamics (RANDALL, 1964). 

The objectives of the present study were to quantify the 
degree to which survival, growth, and fecundity are habitat 
dependent and to determine whether low population den- 
sities in wave-sheltered habitats are due to low recruitment 
or low post-recruit survival. 


MATERIALS anD METHODS 


Field work was conducted between August 1984 and Jan- 
uary 1985 in the Exuma Cays Land and Sea Park, between 
Soldier Cay and Waderick Wells Cay, Bahamas (Figure 
1). The park spans a total of 35 km of the Exuma island 
chain and is 13 km wide. Much of the intertidal zone is 
rocky and offers potential habitat for Cittarum. With winds 
predominantly from the southeast, the eastern shores of 
the islands are exposed to ocean swell from the Exuma 
Sound, while the western shores remain sheltered. 

Four categories of rocky habitat were distinguished along 
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Figure 1 


Map of the study area in the Exuma Cays Land and Sea Park, Bahamas, showing the size-frequency collection 
sites (a-n) and the tag-recapture sites (1-14). Estimates of population density were obtained at all numbered sites 
except 4 and 13. Habitat types are EXP (O), EXSS (A), INT (4), and SH (@). 


this gradient of wave exposure. Exposed (EXP) habitat 
generally showed a 30° slope through the splash zone and 
a nearly vertical drop of 1-1.5 m from the edge of the low 
tide bench (Figure 2). One-meter-high waves were com- 
mon. The substrate in this habitat was highly rugose and 
pitted. The sea urchin Echinometra lucunter and the alga 
Padina (Padina sp.) were abundant on the low tide bench. 
In the upper splash zone, the limpet Acmaea leucopleura 
was abundant. EXP habitat was predominant on eastern 
shores. 

Exposed shallow-sloping (EXSS) habitat contrasted with 
EXP habitat in having a shallow-sloping seabed below the 
low tide level. EXSS habitat was rare and was the only 
other habitat type distinguished on wave-exposed shores. 
Within the study area, EXSS habitat was found on Soldier 
Cay and Halls Pond Cay. 

Intermediate (INT) habitat was found in channels be- 
tween the islands where the effects of ocean swell were 
still present and wave heights of up to 0.5 m were common. 
The habitat had no distinguishing topography, with in- 
tertidal slopes varying widely. Echinometra lucunter was 
less abundant than in EXP habitat. 

Sheltered (SH) habitat was found on the shallow lee- 
ward side of the islands where waves were usually less 


than 0.2 m. Typical were the presence of an undercut in 
the intertidal zone (Figure 2) and less rugosity than in 
EXP habitat. The undercut caused the intertidal zone of 
this habitat type to be shaded for large parts of the day. 
Echinometra lucunter was rare, and the lower-intertidal 
zone was covered by a spongy, sand-laden algal mat. ‘The 
isopod Ligza sp. descended in great numbers into the in- 
tertidal zone at low tide. 

Densities per unit area were obtained for Cittarium and 
predatory gastropods of the genera Thais and Purpura 
(thaids) at four EXP, six INT, and six SH sites. At each 
site, 10 m? of habitat was sampled over the zone of highest 
snail density, following HUGHEs (1971). ANOVA F-tests 
were utilized to compare snail densities between habitat 
types (KLEINBAUM & KUPPER, 1978). All statistical tests 
presented in this report are two-tailed tests. Unless oth- 
erwise stated, statistical comparisons of parameters be- 
tween habitat types were not done by pooling recoveries 
between sites but by considering each site within a habitat 
class as a single replicate. 

Population size-frequency collections were made during 
the peak of the reproductive season (DEBROT, 1990) at 14 
sites in August 1984. Of these, six were EXP, three were 
INT, and five were SH sites. Sample sites were chosen to 
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be at least 150 m apart unless they were separated by a 
sandy beach that prevented dispersal of postlarval indi- 
viduals of this snail. Distribution of the snails was confined 
to a narrow, 1-2 m wide intertidal zone, which greatly 
facilitated sampling. At each site all snails were collected, 
by snorkeling and wading, from as many meters of coast 
as was necessary to fill a standard net bag. 

To quantify survival and growth, a transplant experi- 
ment involving 14 sites was done. A total of 600 small 
snails with a mean shell width of 22.8 + 7.2 mm (95% 
confidence interval) were collected from the ocean side of 
Soldier Cay. Shell width was defined as the widest di- 
ameter across the base of the shell. Oval 2 x 5 mm plastic 
tags were attached to the dry shell surfaces of the snails 
using Marine epoxy cement (Underwater Patching Com- 
pound, Pettit Paint Co., Spring Valley, CA). One hundred 
of the snails were double-tagged to estimate tag shedding 
by the method of GULLAND (1963). The snails were set 
out at four EXP, one EXSS, four INT, and five SH sites. 
Forty snails were released at each site except site 4 (the 
only EXSS site) where 80 snails were released. In addition, 
a total of 101 large (80.0 + 13.8 mm) animals were re- 
leased at three EXP and three SH sites. After approxi- 
mately 168-174 days at liberty, all survivors that could be 
found within a 20-m radius from the point of release at 
each site were collected. 

Handling mortality was estimated by recording the live 
tagged snails seen at each site three days after release and 
assuming that all handling mortality takes effect within 
the first three days. After a period of time, Aż, the fraction 
of snails recovered from the total number of snails released 
can be expressed as RT, = RS, H, where RS, is the recovery 
rate after period Af, assuming no handling mortality, and 
H is handling survival. If natural mortality during the 
three-day period is assumed negligible, then handling sur- 
vival can be calculated as: 


ee — ely RS oe (1) 


where RS,_; is the fraction recovered from snails seen alive 
three days after release. Estimates were tested for differ- 
ences with 100% handling survival using a x? goodness- 
of-fit test (SOKAL & ROHLF, 1981). 

Recovery efficiency was defined as the percentage of 
tagged snails recovered from the site after repetitive sam- 
pling. At each site, tagged snails were removed during low 
tide events until a collection yielded no recoveries. Recovery 
efficiency was calculated separately for small and large 
snails using the DeLury method (RICKER, 1975) on data 
pooled between SH, INT, and EXP sites. 

After a period at liberty, expected tagged-snail popu- 
lation densities should be highest near their point of release 
and should decrease with increasing distance from the 
release point, assuming random movements. Emigration 
was, therefore, quantified by fitting a curve to the number 
of recaptures as a function of distance from the point of 
release (POOLE, 1974). Of the several curves tried, log- 
normal curves of the form Y = a — B In(X) where Y is 


Figure 2 


Vertical shore profiles for (a) EXP and (b) SH habitat types. 
LWS = low tide level, HWS = high tide level (spring tides). 


the number of recaptures at distance X (m) each side from 
the point of release, provided the best fits to the dispersal 
data. Data for habitats with recoveries at multiple sites 
were pooled and weighted for equal representation of each 
site. 

Survival for the period Aż for snails of size 7 was then 
calculated as: 


Sip = Ry /(NioTaF EH) (2) 


where R; is the number of snails of size 7 that were re- 
captured at the end of the period At, No is the number of 
snails of size 7 that were released, 7a is the fraction of tags 
still attached, F; is the fraction of snails of size ¿ inside the 
study site that were also found, Æ is the fraction of snails 
still inside the study site, and H is handling survival. Losses 
due to fishing mortality were considered negligible because 
of the remoteness of the area and the abundance of other 
attractive fishing sites in the area. Survival estimates were 
expressed in terms of annual percent survival (RICKER, 
1975). Data were non-normally distributed so that these 
estimates were compared between habitats using the Krus- 
kal-Wallis test (SOKAL & ROHLF, 1981). 

Dead shells (either empty or occupied by a hermit crab) 
were collected to assess the relative importance of different 
agents of natural mortality. To avoid double-counting, only 
relatively fresh shells with shiny nacre and intact spire 


Page 366 


Table 1 


Population densities of Czttarium and its molluscan pred- 
ators (thaids) in three habitat types: exposed (EXP), in- 
termediate (INT), and sheltered (SH). 


Density (numbers/m?’) 


Nes C. pica Thaids 
Habitat of sites Aver- Aver- 
type sampled age Range age Range 
EXP 4 6.0 3.6-8.5 0.7 0.3-1.1 
INT 6 2.3 0.2-5.9 0.0 — 
SH 6 0.1 0.0-0.5 0.0 — 


tips were used. Shells with extensive breakage were as- 
sumed to have died from crushing predation. Shells with 
either lethal drilling or both lethal drilling and lethal shell 
breakage were assumed to have died due to drilling. Com- 
parisons of the relative importance of various causes of 
natural mortality between habitats were done using the x? 
goodness-of-fit test. The Fisher-Irwin test (x? test-equiv- 
alent for cases with low expected cell-frequencies, HODGES 
& LEHMANN, 1970) was used for equivalent comparisons 
for specific 10-mm size classes. 

Mean shell-width growth increments for the tag-recap- 
ture period were expressed as annual shell growth rates 
and were compared between habitats using an ANOVA 
F-test. Regression analysis was used to estimate the von 
Bertalanffy growth parameters, Læ and K, at sites where 
both large and small snails were released (RICKER, 1975). 
The GULLAND & HOLT method (1959) was used to adjust 
for differences in release time. Confidence intervals were 
obtained using the routine presented by ABRAHAMSON 
(1971). 

Visceral coils of about 50 snails from each sampling site 
were preserved in 70% isopropyl alcohol for determination 
of a gonad index. All snails showing gonadal tissue as 
determined by visual inspection were considered sexually 
mature (RANDALL, 1964). A gonad condition index for 
each specimen was obtained following FEARE (1970). Us- 
ing a color photo of a cross-section of the visceral coil, the 
area corresponding to the gonad was expressed as a fraction 
of the cross-sectional area of the entire visceral coil. Av- 
erage size at first sexual maturity and size-specific fe- 
cundity were compared between habitats using ANOVA 
F-tests. 


RESULTS 


Cittarium population densities were significantly correlated 
with exposure (F413, = 15.85, P < 0.001). Average den- 
sities were highest at EXP sites (6.0 snails/m?) and lowest 
at SH sites (0.1 snails/m?) (Table 1). Predatory thaid 
gastropod densities were similarly correlated with expo- 
sure (F913, = 19.84, P < 0.001). Thaids were abundant 
at EXP sites (0.7 thaids/m?) but virtually absent at the 
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Table 2 


Cittarium releases, percent recovery, and percent handling 
survival (H) for small (* = 22.8 mm) and large (x = 80.0 
mm) snails. 


Seen 3 days 


Total releases after release 


Num- Percent Num- Percent 


ber re- recov- ber recov- 

leased ered seen ered* H 
Small snails 596 16.4 337) 17.5 94.0 
Large snails 101 23.8 63 27.0 88.0 


* Percent recovered from those snails that were observed to be 


alive 3 days after release (RS,; in equation 1). 


INT and SH sites examined (Table 1). At EXP sites, 
Thais deltoidea constituted 64% of the 96 thaids seen, while 
T. rustica and Purpura patula constituted 24% and 12%, 
respectively. Sixteen of the 18 thaids (89%) recorded over 
the course of more than two years at INT and SH sites 
were T. deltoidea, while two (11%) were P. patula. 

Few population modes were evident in the size-fre- 
quency data obtained (Figure 3). Population structures at 
EXP sites were dominated by small snails (z.e., 15-30 mm 
shell width). Despite large sample sizes (n = 174-446), 
low abundances of large snails precluded meaningful dis- 
tinction of more than one population mode. Lack of snails 
in the larger size classes suggested a high rate of extinction 
of population modes due to low survival or low growth 
rates at EXP sites. In contrast, population structures at 
SH sites were dominated by large snails near the maximum 
size attained by the species. Indistinct modes at SH sites, 
due to small sample sizes (n = 56-118) and combined with 
the fact that population modes near Le tend to merge 
(SAINSBURY, 1982), precluded more detailed analysis of the 
size-frequency data for survival and growth estimation. 
The size-frequency data, nevertheless, suggested important 
differences between habitats in growth or survival. 

A total of 19 double-tagged shells were recovered with 
tags still attached. Of these, 9 had one tag and 10 had two 
tags. Annual tag shedding was estimated at 55% and did 
not differ significantly from an estimate obtained in a 
concurrent study (38%; n = 43) using the same tags (DE- 
BROT, 1990). Therefore, the combined estimate of 43%/ 
year was used for the purpose of survival estimation. 

Handling survival was not significantly less than 100% 
for either the large (x? = 0.22, d.f. = 1, P > 0.5) or small 
snails (x? = 0.18, d.f. = 1, P > 0.5) released (Table 2). 
Recovery efficiency for small snails did not differ signifi- 
cantly between habitat types (x? = 2.40, d.f. = 3, P > 
0.50). By the end of recovery, 96% of the small tagged 
snails present within 20 m from the points of release at 
INT and SH sites had been recovered (Figure 4). Recovery 
efficiency for large snails was estimated at 97%. For both 
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Figure 3 


Size-frequency distributions for Cittarium collected at the 14 size- 
frequency collection sites (a-n), grouped by habitat type. 


small and large snails, a 100% recovery efficiency was 
assumed. 

Under ideal conditions (7.e., large sample size, no mi- 
crohabitat selection, and uniform behavior by all members 
of the population), snail recaptures should decrease with 
increasing distance from their point of release. In general, 
recaptures decreased with increasing distance from their 
point of release (Figure 5). At the EXSS site an estimated 
94% of the snails originally released had remained within 
a 20-m radius from the point of release. At INT sites 67% 
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Figure 4 


Cumulative number of tagged Cittarium recaptured 168-174 days 
after release versus number of successive collections, for all tag- 
recapture sites combined. Also shown are the asymptotically es- 
timated numbers of large (----- ) and small (~ ~- -) snails present 
within 20 m from their points of release. 


of the snails had remained within a 20-m radius while at 
SH sites 85% had. After scaling for differences in total 
number of recoveries per site, habitat differences in dis- 
persal were significant (x? = 15.41, d.f. = 4, P < 0.005). 

No dispersal correction factor was calculated for EXP 
habitat because of few recoveries. However, the two re- 
coveries of small shells made at EXP sites were within 1 
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Figure 5 


Percentage recoveries of live tagged Cittarium as a function of 
the distance from their points of release, after 168-174 days at 
large. Fitted curves are of the form Y = a — B In(X) where Y is 
the number of recaptures and X is the distance from the point 
of release. 
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Table 3 


Tagged-shell releases and recoveries for small Cittarxum (x = 22.8 mm) by site, and the corresponding growth and survival 
estimates. Dashes indicate data not available. Habitat types: exposed (EXP), exposed shallow sloping (EXSS), intermediate 
(INT), and sheltered (SH). 


Number of snails 


Habitat type Site number Released Recovered 

EXP 1 40 2 
2 40 0 

3 40 0 

5 40 0 

EXSS 4 80 22 
INT 6 40 6 
7 40 0 

8 40 2 

9 40 13 

SH 10 40 8 
11 40 4 

12 36 10 

13 40 19 

14 40 12 


Survival rate Growth rate 


Percent recovered (%/year) (mm/year) 

5 0.7 7.9 

0 <0:2* — 

0 =<02% — 

0 <02 — 
Average: 0.2 TO) 
Biles 12.7 10.2 
15.0 YD) 15.2 

0.0 <0.2* — 
5.0 0.7 15.2 
32:5 375 : 25.0 
Average: 11.3 18.4 
20.0 8.0 22.6 
10.0 1.8 19.4 
27.8 16.1 22.4 
47.5 50.7 31.2 
30.0 19.0 16.9 
Average: 19.1 22.4 


* Due to the small number of releases and the fact that fractional recoveries are not possible, the estimates of survival are not 
continuously distributed. For the cases with zero recoveries, mortality estimates are presented as less than (<) the value that would 
have been obtained with one recovery. In all such cases, a zero survival value was used in calculating average annual survival for each 


habitat type. 


m from the point of release. The low number of recoveries 
at EXP sites is not believed to be due to unusually high 
levels of emigration because searches far beyond 20 m from 
the point of release yielded no additional recoveries. Sur- 
vival estimates for EXP sites were, therefore, derived using 
the dispersal correction factor for INT sites, the most 
conservative choice. 

Survival estimates for small snails differed significantly 
between the habitats for which multiple estimates were 
obtained (P = 0.032, Kruskal-Wallis test). Lowest average 
survival was at EXP sites (0.172%/year) where only 2 of 
the 160 snails released were again recovered. Average sur- 
vival for small snails was much higher at other sites, av- 
eraging 12.7%/year, 11.3%/year, and 19.1%/year for 
EXSS, INT, and SH sites, respectively (Table 3). For the 
large snails released, average survival estimates were 6.2% / 
year at EXP sites and 29.8%/year at SH sites (Table 4) 
but, with a low number of replicates per habitat type, did 
not differ significantly (P > 0.10, Mann-Whitney U test). 

Few large dead shells were collected from EXP sites 
(Figure 6), mirroring their low abundance in the live snail 
collections. At both EXP and INT sites, most dead shells 
were classified as lethally crushed (50% and 42%, respec- 
tively). At EXP sites, drilled shells constituted the second- 
largest fraction, while undamaged shells were the least 
abundant of the three fractions. In contrast, at INT sites, 
undamaged shells constituted the second-largest fraction 


while drilled shells were the least abundant. Higher frac- 
tions of drilled and crushed shells at EXP sites suggested 
higher predation mortality than at INT sites (x? = 20.10, 
d.f. = 2, P < 0.005). However, due to a tendency for larger 
shells to have a lower incidence of predatory shell damage 
(Table 5) and due to the difference between habitats in 
dead-shell size structure, size-specific comparisons were 
necessary. For snails in the size classes for which com- 
parable samples were obtained (the 30-40 mm and 40- 
60 mm size classes), P-values of less than 0.06 suggested 
higher predation mortality in EXP (Table 5). 

Average growth rates for small snails were 7.9 mm, 10.2 
mm, 18.4 mm, and 22.4 mm/year for EXP, EXSS, INT, 
and SH sites, respectively (Table 3). Lowest growth rates 
were in EXP and EXSS habitats, both of which had high 
snail densities. Small sample sizes precluded statistical 
comparison of more than two habitats unless estimates 
from EXP and EXSS habitats were pooled to form a third 
habitat class. A comparison of growth rates between INT 
and SH habitats showed no significant differences (Fa, 
= 1.032, P = 0.349), whereas including a third habitat 
class (with pooled EXP and EXSS estimates) did result 
in significant habitat differences (Fy 7) = 4.918, P = 0.046). 
For the more limited release of large tagged snails, average 
growth rates appeared uniformly low among habitat types 
(Table 4). 

Von Bertalanffy growth parameter estimates for EXP 


A. O. Debrot, 1990 
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Table 4 


Tagged-shell releases and recoveries for large Cittartum (¥ = 80.0 mm) by site, and the corresponding growth 
and survival estimates. 


Number of snails 


Survival rate Growth rate 


Habitat type Site number Released Recovered Percent recovered (%/year) (mm/year) 
EXP 2 17 3 17.6 10.2 4.5 
3 17 2 11.8 4.3 2.6 
5 17 2 11.8 4.3 8.1 
Average: 6.2 Sl 
SH 10 17 6 35.3 26.9 3.4 
13 16 8 50.0 56.5 1.9 
14 17 3 17.6 6.1 2.1 
Average: 29.8 2.6 


sites also suggested low growth rates compared to SH sites. 
Using midpoint parameter estimates (Table 6), the age of 
a 65-mm snail from SH habitat would be about four years, 
while that of a 65-mm snail from EXP habitat would be 
about nine years. 

Average shell width at first maturity was smallest at 
EXP sites (¢ = 16.7 mm; range, 14.5-17.8), intermediate 
at INT sites (x = 20.2 mm; range, 18.4-21.4), and largest 
at SH sites ( = 24.2 mm; range, 22.8-27.9), and differed 
significantly between habitats (Fy 49) = 23.37, P = 0.002). 
In addition, fecundity levels were consistently lowest at 


EXP sites (Table 7). 


DISCUSSION 


The tagging experiment indicated highest rates of mor- 
tality for Cittartum in EXP habitat. The higher density of 
molluscan predators and the higher proportion of lethal 


predatory shell damage, as compared to less exposed hab- 
itat, may account for these observations. Growth rates, size 
at first sexual maturity, and fecundity were the least in 
EXP habitat. Higher Cittartum density may reduce per 
capita food availability in EXP habitat. However, other 
factors, such as differences in wave stress, air and water 
temperatures, and floral food value, may also be important 
factors affecting growth rate, size of maturation, and fe- 
cundity. 

In contrast, EXSS habitat had high densities of large 
animals. Growth and survival rates were measured at one 
EXSS site (site 4). Growth, based on 22 recaptures, was 
similar to that at other high-density EXP sites and was 
slow compared to growth rates at SH and INT sites. The 
survival estimate at the EXSS site was much higher than 
at EXP sites, suggesting that shallow topography may 
reduce mortality rates. In particular, dislodgement by waves, 
which is known to subject snails of wave-exposed habitat 


Table 5 


Number of Cittartum shells in each of three shell-damage categories for shells collected from EXP and INT habitats. 


Habitat Shell damage 

type category 0-20 21-30 31-40 

EXP Drilled 26 61 41 
Crushed 14 109 43 
Undamaged 10 24 9 
Totals 50 194 93 

INT Drilled 1 2 3 
Crushed 0 10 7 
Undamaged 0 1 11 
Totals 1 13 21 
BS — 0.952 0.004 


‘ Percent of 
Shell size (mm) aa 
41-50 51-60 >61 All sizes servations 
11 3 1 143 37 
15 11 1 193 50 
3 2 5 56 13 
29 16 7 389 100 
4 3 9 22 25 
11 4 T5 37 42 
5 5 y 29 33 
20 12 21 88 100 
Coa a 
0.059 0.194 <0.005 


* P-values compare size-specific importance of predatory mortality (combined crushed and drilled) and nonpredatory mortality 
(undamaged) between habitats. Fisher-Irwin tests except the comparison for all shell sizes combined. The latter was done using a x? 


goodness-of-fit test. Dashes indicate test not made. 


